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Relativistic energy density functionals have become a standard framework for nuclear 
structure studies of ground-state properties and collective excitations over the entire 
nuclide chart. We review recent developments in modeling nuclear weak-interaction pro¬ 
cesses: charge-exchange excitations and the role of isoscalar proton-neutron pairing, 
charged-current neutrino-nucleus reactions relevant for supernova evolution and neu¬ 
trino detectors, and calculation of /3-decay rates for r-process nucleosynthesis. 
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1. Introduction 

Nuclear weak-interaction processes play a crucial role during various phases of su¬ 
pernova evolution and in the associated production of chemical elements.®^ In the 
presupernova phase, electron capture and /3-decays produce neutrinos and change 
the number of electrons in a way that directly influences the collapse dynamics and 
subsequent explosionP The r-process nucleosynthesis, responsible for the synthe¬ 
sis of about half the elements heavier than iron, is governed by neutron capture 
and ,5-decays. Neutrino-induced reactions on nuclei also play an important role in 
supernova environment. Neutral-current neutrino-nucleus scattering in stellar envi¬ 
ronment determines the rate of cooling by neutrino transport.®^ Inelastic neutrino- 
nucleus scattering has been included in supernova simulations as a novel mode of 
energy exchange between neutrinos and matter Furthermore, charged-current 
neutrino-nucleus reactions are important in the environment of exploding massive 
stars, considered as one of the possible sites for r-process nucleosynthesis. Neutrino- 
induced reactions in explosive supernova nucleosynthesis significantly contribute to 
the production of not only radioactive, but also long-lived nuclides.^^ 


1 
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As emphasized in Ref. |3l a real breakthrough in our understanding of the in¬ 
ner workings of supernova explosions, based on self-consistent models with all the 
relevant microphysics included, has not been achieved yet. The necessary nuclear 
input is obtained either from a variety of independent sources based on different ef¬ 
fective interactions and approximations, or it is based on relatively crude estimates 
that cannot provide a realistic description of key weak-interaction processes.!^ Self- 
consistent modeling of heavy-element nucleosynthesis presents a particularly diffi¬ 
cult challenge, both because of the nature of the r-process and the large amount 
of nuclear data input involved. As the detailed description of this process requires 
a knowledge of the structure and decay properties of thousands of nuclei, most of 
which are outside the reach of present-day experimental facilities, it is necessary to 
develop a universal self-consistent theoretical framework that can systematically be 
used over the entire nuclear chart. One of the key phenomena involved is /3-decay, 
and the closely related /3-delayed neutron emission, as it determines the rate of the 
r-process and directly influences the final elemental abundances. We note, however, 
that there are few systematic studies of the impact of first-forbidden transitions on 
the total decay rates of neutron-rich nuclei. 

In this article we present an overview of current developments in self-consistent 
modeling of weak-interaction processes and associated nuclear transitions based on 
the framework of relativistic nuclear energy density functionals (RNEDF). The fo¬ 
cus is on charged-current neutrino-nucleus reactions and /3-decays. Many results 
of alternative theoretical approaches to ,5-decay rates, neutrino-nucleus reactions 
and their role in supernova dynamics and nucleosynthesis have been reviewed in 
Refs.fTlElISl ITniT^ and references therein. Over the last decade significant progress 
has been achieved in modeling ground-state nuclear properties, excitations and as¬ 
sociated weak-interaction phenomena using the RNEDF framework. Recent studies 
include the analysis of collective nuclear phenomena, including giant resonances and 
exotic modes of excitation, E^HUl their application in constraining the neutron 

skin of neutron-rich nuclei,I^^J^ symmetry energy and neutron star properties.!^ 
RNEDF-based models have also been used in calculations of transition matrix ele¬ 
ments that contribute to a variety of weak-interaction processes, including ,5-decay 
of r-process nuclei,!22l[^ muon captur^^ and stellar electron capture,!^ neutral- 
currentP^ and charged-current neutrino-nucleus reactions.!2^1!^ The investigation 
of charge-exchange transitions in nuclei and comparison with available data pro¬ 
vide important benchmark tests for microscopic models of weak-interaction pro¬ 
cesses 

The paper is organized as follows. The theory framework is briefly outlined in 
Sec. [21 In Sec. [3] we present a benchmark test for charge-exchange excitations in 
®®Fe. Sec. [4] explores the effect of isoscalar pairing interactions on Gamow-Teller 
transitions in light and medium-mass nuclei. In Sec. [5] charged-current neutrino- 
nucleus cross sections are considered, and in Sec. [6] results of large-scale calculations 
are analyzed. /3-decay half-lives and results of systematic calculations are discussed 
in Section [71 Sec. [8] includes a brief summary and concluding remarks. 
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2. Theory framework 

The analysis of collective excitations and nuclear weak-interaction processes is based 
on the framework of relativistic nuclear energy density functionals (RNEDFs). 
At the mean-field level ground-state calculations of open-shell nuclei are imple¬ 
mented using the relativistic Hartree-Bogoliubov (RHB) model, ^ based on density- 
dependent meson-nucleon effective interactions!^ A detailed descripton of the for¬ 
malism and the corresponding computer codes are included in RefP^ In this work 
we employ the relativistic functional DD-ME^^ in the particle-hole channel, and 
pairing correlations in open-shell nuclei are described by the finite-range Gogny 
interaction DISp^ In the small-amplitude limit collective nuclear excitations rele¬ 
vant for weak-interaction processes are calculated using the relativistic quasiparticle 
random-phase approximation (RQRPA) The RQRPA is formulated in the canon¬ 
ical single-nucleon basis of the RHB model, and the residual interaction is derived 
from the same relativistic functional used to compute ground state properties. An 
important advantage of this framework is that it employs in a consistent way the 
same effective interaction both in the RHB equations that determine the canonical 
basis, and in the matrix equations of the RQRPA. 

3. Charge-exchange excitations 

Gamow-Teller(GT^) excitations in nuclei correspond to spin-flip and isospin-flip 
transitions, involving both the spin and isospin operators (7t±. Recently GT^ tran¬ 
sitions in 54’^®Fe and ^®>™Ni have been explored in a comparative study based 
on the RNEDFs, Skyrme energy density functionals, and the shell model.^ As a 
benchmark case. Fig. [T] displays the GT“ strength distributions for ®®Fe, computed 
using two different theory frameworks: the relativistic QRPA based on the func¬ 
tional DD-ME2 and the shell model (GXPFIJ). The data from {p,n) reaction^^ 
are also included for comparison. The calculated transition strengths are folded by a 
Lorentzian function of width F=0.5 MeV. More details about the shell model Hamil¬ 
tonian GXPFIJ and its implementation in the pf shell nuclei are given in Refs.^^^®^ 
Although the QRPA model calculation includes only 2qp configurations and, there¬ 
fore, cannot reproduce the detailed structure of excitation spectra obtained by the 
shell model, one finds a reasonable agreement with the global features of the GT“ 
transition strength distribution. We note, however, that the measured overall GT^ 
transition strengths in and can only be reproduced by quenching the 

free-nucleon axial-vector coupling constant g^.l^This is true for all models and the 
respective parameterizations employed in the study: i) RNEDF (DD-ME2), RPA 
based on Skyrme functionals (SGH, SLy5) and iii) the shell model (GXPFIJ). 
The quasiparticle RPA calculations include the quenching of the free-nucleon axial- 
vector coupling constant gA = 1.262 g a = 1, corresponding to the quenching 

factor 0.8 in the GT transition operator. A quenching factor 0.74, used in shell 
model calculations, is adapted to the effective interaction and model space under 
consideration.!^ In the case of ®®Fe, the RNEDF QRPA calculation with = 1 
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yields the B(GT“) within experimental uncertainties!^ The quenching of gA, how¬ 
ever, should be considered with caution because of the well known problem of 
missing GT strength, either due to excitations that involve complex configurations 
at higher excitation energies,^ or excitations that include non-nucleonic degrees of 
freedomP^ 



Fig. 1. The Gamow-Teller (GT“) transition strength distribution for ®®Fe, shown as a function 
of excitation energy in the final nucleus. The RQRPA results based on the RNEDF DD-ME2 are 
compared to the shell model calculations (GXPFIJ),!^ and available data from (p, n) reactions.!^ 


4. Gamow-Teller transitions and the isoscalar pairing interaction 

Microscopic descriptions of open-shell nuclei based on energy density function¬ 
als necessitate the inclusion of an effective pairing interaction. The isovector 
(T = 1,5 = 0) part of of this interaction accounts for the odd-even staggering 
in the separation energies, and for the gap in the excitation spectrum of even-even 
nuclei.123®^ In the relativistic framework, standard mean-field methods that include 
pairing correlations in the description of ground-state properties of open-shell nu¬ 
clei are the relativistic Hartree -I- Bardeen-Gooper-Schrieffer (BGS) model,!^ and 
relativistic Hartree-Bogoliubov (RHB) model.^ A fully self-consistent calculation 
of excitations in open-shell nuclei involves the inclusion of particle-particle corre¬ 
lations not only at the ground-state level, but also in the residual interaction of 
the quasiparticle RPA equations.^® In the case of charge-exchange transitions, the 
particle-particle interaction includes not only the usual isovector (T = I) part, 
which has the same form as the pairing interaction used for ground-state calculation 
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(e.g., the Gogny pairing), but also the isoscalar T = 0 channel, for which one could 
adopt the following form: 


2 

Vi2 = ns=i_T=o ) (1) 

i=i 

where ns=i_T=o projects onto states with S' = 1 and T = 0. The standard choice 
for the ranges of the two Gaussians is /ri=1.2fm and ^2=0.7 fm, together with 
the relative strengths gi = 1 and g 2 = —2 so that the force is repulsive at small 
distancesP^ The remaining free parameter is Vq, the overall strength that can only 
be constrained by data on nuclear excitations and/or decays. As pointed out in 
several HFB+QRPA calculations, by adjusting the strength parameter of isoscalar 
pairing interaction to selected data, one can improve the description of /3-decay half- 
lives in neutron-rich nuclei relevant for the r-process nucleosynthesis (see Sec. [7] for 
more details).^^^!^^ 

In the study of Ref. |49] it has been shown that the T — 0 particle-particle inter¬ 
action leads to considerable modification of the Gamow-Teller strength distribution 
in Sn isotopes. In particular, the energy spacings between the main GT resonance 
peaks depend on the strength of T = 0 pairing, and the low-energy GT transi¬ 
tion strength becomes more pronounced by increasing Vq. Because of the complex 
dependence of GT transition strength on the T = 0 pairing strength in various 
mass regions, it is difficult to determine a unique global value of Vq from available 
experimental results. Since the modeling of /3-decay rates crucially depends on the 
low-energy part of GT excitation spectra,!^ a possible approach is to employ ex¬ 
perimental /3-decay half-lives to constrain Vq locally for specific groups of nuclei, or 
to introduce a smooth mass dependence for the strength parameter (Sec. [7]). 

The role of T = 0 pairing for GT states in N = Z (A = 48 — 64) nuclei was 
investigated using the Hartree-Fock-Bogoliubov + QRPA model based on Skyrme 
functionals.!^ It has been shown that the low-energy excitation strength is enhanced 
by T = 0 pairing and that N = Z nuclei could provide important information about 
this channel of the effective interaction in medium-heavy and heavy mass nuclei.!^ 
Recently the GT transition strength distributions were explored in {^He, t) charge- 
exchange reactions with f-shell nuclei,!^ resulting in a dominant concentration of 
GT strength in low-energy states, in particular for ^^Ga, ^®Ti, and ®°Gr. As shown 
in Refs. EH [551 the isoscalar spin-triplet particle-particle interaction mixes largely 
the {vjy = Z -I- 1/2 —)• 7r/< = I — 1/2) configurations into the low-energy states, and 
this mixing plays an important role in the formation and evolution of collectivity 
of low-energy GT states. 

In Fig. [2] we plot the RHB-I-RQRPA GT“ transition strength distributions for 
i 8 , 20 , 22 q self-consistent ground states are computed using the RHB model 

(DD-ME2 functional and Gogny pairing in the T = 1 channel), and the transition 
strengths are obtained in the relativistic QRPA that includes both the T = 1 and 
T = 0 pairing interactions. In the T = 0 case, the strength parameter is varied in the 
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Fig. 2 . Gamow-Teller {GT“) strength distributions for 18,20,22 q calculated with the 
RHB+RQRPA model (DD-ME2 functional and Gogny pairing in the T = 1 channel). The strength 
of the T = 0 particle-particle interaction Eq. |[T} varies from Vq = 0 to Vq = 300 MeV. 


range of values from Vb = 0 to Vb = 300 MeV. For the main GT“ state at 10.8 
MeV (Vb=0 MeV) corresponds predominantly to the i'ld ^/2 —t 7rl(i3/2 transition. 
However, by increasing the strength of the T = 0 pairing interaction the structure 
of the main GT“ peak becomes more complex. For example, for Vb=300 MeV 
the main peak is shifted to 7.5 MeV, and the dominant configurations 
'^1<^3/27 7rl(i5/2, v^S\i2 7'’ld3/27 i^l'^3/2 7r2s3y2! v^S\i2 7'’2si/2i and 

v\d^l 2 7rld5/2, contribute with 65%, 10%, 10%, 6%, 5%, and 1%, respectively, to 
the QRPA amplitude X^—Y^. We note that for all the additional configuration pairs 
except the last one, the occupation probabilities of both single-particle states are 
rather small, i.e. the contribution of these particle-like configurations is due to the 
T = 0 pairing interaction in the QRPA. Because of destructive interference between 
the i'ld ^/2 7rld3/2 and particle-like configurations, the overall GT“ strength is 

reduced for larger values of Vb- As shown in Fig. [2] for the oxygen isotopes, the 
GT“ transition strength is sensitive to the neutron excess: from to the 
main peak is shifted to lower energy by 5.4 MeV. While for (Vb = 0 MeV) the 
structure of the main peak is dominated by a single transition, for the QRPA 
amplitude of the main peak at 5.42 MeV is almost equally distributed between two 
configurations: vld5/2 7rld3/2 and uld3f2 7rlci3/2. 

Figure [3] displays the calculated GT“ strength distributions for ^^Ca. Similar 
to the previous case, the dependence of transition spectra on the T = 0 particle- 
particle interaction Eq. m is explored by varying the pairing strength. In the 
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case without isoscalar pairing (Vb = 0), one finds a dominant peak at 16.2 MeV 
that corresponds to the transition v\fi /2 —t 7rl/5/2. Another prominent state is 
calculated in the low-energy region at 8.9 MeV {i'lfr /2 ^^f 7 / 2 )- For Vb = 0 
MeV, the B(GT) values for the main transitions {vj> —>■ 7rj>) and {vj> —>■ 7rj<) 
exhaust 34% and 53% of the GT sum rule 3{N — Z) = 6 , respectively. However, 
recent data from {^He,t) reactions indicate that 80% of the GT strength is carried 
by the lowest GT state.^The T = 0 particle-particle interaction provides a possible 
mechanism to shift the calculated strength toward lower excitation energies.^^ 
This effect is demonstrated in Fig. [3] for the range of values of the T = 0 pairing 
strength parameter Vb = 50 — 300 MeV. By increasing Vb the peaks are lowered in 
energy, the B(GT) value of the low-energy state is considerably enhanced, whereas 
the peak of the high-energy state decreases. For Vb = 300 MeV the low-lying state 
exhausts 59% and the high-lying state 33% of the GT sum rule. The origin of 
the discrepancy with experimental results, which show no evidence of pronounced 
high-energy GT transitions,^ remains an open problem. As pointed out in Ref. 1551 
an additional effect could be due to contributions of the tensor force. However, 
in the case of '*^Ca this effect is negligible in the low-energy states, whereas the 
high-energy states are shifted downward in energy by 1-2 MeV.^ 



Fig. 3. Same as the caption to Fig. m but for '‘^Ca. 
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5. Charged-current neutrino-nucleus reactions 

In charged-current neutrino-nucleus reactions an incoming electron neutrino (ve) 
induces a charge-exchange in the target nucleus X{Z,N), 


zXn 


z+iXj^_i -|- e 


( 2 ) 


The general formalism for calculating neutrino-nucleus cross sections, derived using 
the weak-interaction Hamiltonian in the current-current form, is given in detail in 
Refs. [5S11571 The cross section includes the transition matrix elements of the charge 
Xij, longitudinal £j, transverse electric and transverse magnetic 7^^*^ mul¬ 
tipole operators, between the initial \ Ji) and final \ Jf) nuclear states.!^ A complete 
description of inelastic neutrino-nucleus scattering in the energy range of supernova 
neutrinos requires the inclusion of higher multipoles jW^ Although the contribution 
from higher-order multipole transitions to the cross section is rather small at low 
incoming neutrino energies, these cannot be neglected at neutrino energies of tens 
of MeV. 

When interacting with nuclei, neutrinos induce transitions from the initial state 
of a target nucleus (i.e. the ground state J'^=0+) to a large number of excited states 
that are allowed by conservation of energy and selection rules. The contribution 
of various multipole transitions to the neutrino-nucleus interaction cross sections 
essentially depends on the energy of incoming neutrinos. Figure U displays the 
contributions of the multipole transitions = 0* — 4^ to the inclusive cross 
section for the ®®Fe(r'e, e“)®®Co reaction, at selected neutrino energies = 40, 
60 and 80 MeV. In the RHB-I-RQRPA (DD-ME2) model calculation the value 
of axial-vector coupling constant is used for all multipole operators. The 

partial contributions are compared with those of hybrid calculation that employs the 
shell model (GXPFIJ) for the I’*' channel, whereas the RPA based on the Skyrme 
functional SGII is used for higher multipoles.^ For lower neutrino energies, Ei, < 
40 MeV, the dominant contribution to the cross section originates from J'^=l+ 
transitions. However, with increasing neutrino energy higher multipole transitions 
start contributing to the cross section, as shown for E^=QQ MeV and 80 MeV. We 
note that the two approaches predict a similar structure and distribution of the 
relative contributions from various multipoles. 

Experimental results for inelastic neutrino-nucleus cross sections are rather lim¬ 
ited. For two target nuclei ^^G and ®®Fe, data are available from the LAMPF,!^ 
KARMEISf^ and LSNEI^^HsII collaborations. To compare model calculations with 
the data, cross sections are averaged over fluxes which depend on the neutrino 
source. In the case of t'e, the Michel flux from the decay at rest (DAR) of muons ) 
is usedP^ f{Eve) = — 2E,^J. The flux averaged cross section for 

the target has been calculated using nonrelativisticP^^^ and relativisticP^^^I^ 
energy density functionals, the projected quasiparticle random phase approxima¬ 
tion,!^ and the shell model.lS^The cross sections for the reaction ®®Fe(r'e, e )^®Co 
were recently analyzed to estimate the systematic uncertainty in modeling neutrino- 
nucleus cross sections for medium mass nuclei.!^ The results were compared with 
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Fig. 4. Contributions of the multipole transitions J'^ = 0^ —4^ to the inclusive cross section for the 
^^Fe{ue, e“)^®Co reaction, at E,y^ = 40, 60 and 80 MeV. The results correspond to RHB+RQRPA 
(DD-ME2) calculations, and the shell model (GXPFIJ) (for the 1'*" transition) plus the RPA 
(SGII) for higher multipoles.^^ 


the experimental uncertainty of the the KARMEN data. The calculated i/g-^^Fe 
cross sections, averaged over the Michel spectrum, are listed in Table 1. In addition 
to the RHB+RQRPA model calculation based on the density functional DD-ME2, 
results are shown for the shell model (GXPFIJ) (for 1+ transitions) plus the RPA 
(SGII) for higher multipoles,l^l^ RPA based on the Landau-Migdal force, ^ the 
QRPA(SIII)P^and QRPA based on the G-matrix formalism.l^By using these meth¬ 
ods, the cross section mean value and uncertainty are obtained: < a >t/i=(258±57) 
xio-42cm2, in excellent agreement with the data from the KARMEN collabora¬ 
tion: < (7 >exp=(256±108±43) X 10“"'^^cm^. At present, the systematic theoretical 
uncertainty of the calculated cross section appears to be smaller than the corre¬ 
sponding experimental value. 

In the environment of high neutrino fluxes that occur in core-collapse supernovae 
or neutron star mergers, neutrino-nucleus reactions are particularly important for 
the process of nucleosynthesis. In addition, neutrino detectors are based on inelastic 
neutrino-nucleus scattering and, to provide reliable predictions of supernova neu¬ 
trino induced events, it is essential to be able to compute neutrino-nucleus cross 
sections on a quantitative level. The supernova neutrino fluxes can be obtained 
from core-collapse supernova simulations.^ However, for the purpose of testing 
the sensitivity of the underlying models of nuclear structure and neutrino-nucleus 
interactions, a simplified neutrino flux is employed, described by the Fermi-Dirac 













May 29, 2015 0:48 WSPC/INSTRUCTION FILE collexc'DVMay21 


10 N. Paar, T. Marketin, D. Vale, and D. Vretenar 


Table 1. Comparison of the inclusive t'e-®®Fe cross sections 
averaged with the Michel flux. 


< fT > [10-'‘2cm2] 

RNEDF (DD-ME2) 

263 

SM(GXPFIJ) + RPAfSGII^EH 

259 

RPA (Landau-Migdal)^ 

240 

QRPA (siii;@3 

352 

QRPA(G-matrix;|lIl 

173.5 

Theoretical average 

258±57 

Exp. (karmen;^ 

256±108±43 


spectrum, 

/(Ejy) = nrT ’ 

exp [{E„/T) - aj + 1 

where T denotes the temperature, and a is the chemical potential. As shown in 
Ref. 130! for the set of inclusive neutrino-nucleus cross sections calculated using the 
RNEDF framework, the implementation of more realistic fluxes from core-collapse 
supernova simulations is straightforward. 

For a quantitative analysis of total cross sections, and also partial cross sec¬ 
tions for neutrino-induced particle knockout, in a first step the neutrino-induced 
excitation spectrum in the daughter nucleus has to be computed. In Ref. [57| it has 
been shown that this spectrum provides a basis for the implementation of statis¬ 
tical model codes (e.g. SMOKER), to determine for each final state the branching 
ratios into various decay channels. In Fig. [5] we plot the dominant multipole con¬ 
tributions to the RHB-hRQRPA (DD-ME2) excitation spectrum induced in the re¬ 
action ^®Fe(z/e, e“)^®Co. The cross sections that correspond to the transition from 
the ground state of the target nucleus to each excited state in the daughter nu¬ 
cleus are averaged over the supernova neutrino flux Eq. (|3|), for T = 2,4, and 6 
MeV, and a = 0. Fig. [Sjshows that, by increasing the temperature that determines 
the neutrino flux, the distribution of the neutrino-induced excitation spectrum is 
shifted toward higher energies. In addition, at higher temperatures the structure 
of the excitation spectra becomes rather complex due to neutrinos in the extended 
tail of the Fermi-Dirac distribution. While for T = 2 MeV only a few 1+ and O'*" 
states contribute to the reaction cross section, at T = 6 MeV a large number of 
states of various multipoles yield non-negligible contributions. Figure [6| displays 
the same analysis, but for the case of inelastic antineutrino - ®®Fe scattering. For 
temperatures in the range T = 2 — 4 MeV a single low-lying 1+ state gives the 
main contribution to the cross section, whereas at higher temperatures significant 
contributions arise from excited states at higher energies. 

The flux-averaged inclusive cross sections for the reaction ®®Fe(z^e, e“)®®Co, eval¬ 
uated at different temperatures in the interval T = 2 — 10 MeV, a = 0, are displayed 
in Fig. [71 The results of RHB-fRQRPA (DD-ME2) calculations are compared to 
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Fig. 5. Neutrino-induced excitation spectrum for the reaction ^®Fe(z>'e 5 e“)^®Co, averaged over the 
Fermi-Dirac distribution Eq. ((S]) with T = 2, 4, and 6 MeV, q = 0. The RHB+RQRPA calculations 
are based on the EDF DD-ME2. The excitation energies are given with respect to the ground state 
of the target nucleus. 


those obtained with the shell model + RPA (GXPFIJ + SGII),l^the hybrid model 
of Ref.,and estimates from Ref. |3] (Woosley et ah). As shown in the figure, an ex¬ 
cellent agreement is obtained for the RHB-I-RQRPA, shell model and hybrid model, 
with small deviations noticeable only at the high-end temperature T = 10 MeV. 
At temperatures above T = 4 MeV the results of Ref. H] display deviations with 
respect to more recent microscopic calculations. 


6. Large-scale calculations of neutrino-nucleus reactions 

Inelastic neutrino-nucleus scattering plays an important role in supernova evolu¬ 
tion and nucleosynthesis.^ During the star collapse phase, charged-current (i/e,e“) 
reactions on nuclei are suppressed because of Pauli blocking of the final electron 
phase space.^ However, in the environment of exploding massive stars or neutron 
star mergers, considered as possible sites for the r-process nucleosynthesis, charged- 
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Fig. 6. Same as the caption to Fig. [5] but for the antineutrino induced reactions on ®®Fe. 


current neutrino-nucleus reactions could contribute to the production of heavier 
isotopes. Since a complete modeling of neutrino-induced reactions requires the in¬ 
clusion of not only Gamow-Teller type of transitions but also contributions from 
higher multipoles, including a large pool of target nuclei can be computationally 
very demanding. Even though the shell model provides an accurate description of 
ground state correlations and low-lying excited states, the calculation of high-lying 
excitations and/or higher multipole can be a formidable task, making this approach 
applicable essentially only to allowed transitions in light and medium-mass nuclei. 
For systematic studies of neutrino-nucleus cross sections over the entire nuclide 
chart, microscopic calculations must therefore be implemented using a framework 
based on universal energy density functionals. 

A consistent set of models based on the RNEDF framework was recently em¬ 
ployed in large-scale calculations of charged-current neutrino-nucleus cross sections, 
for a set of even-even target nuclei from oxygen to lead (Z = 8 — 82), and with the 
neutron number spanning the interval N = 8 — 182.^^ These extensive calculations 
include allowed and forbidden transitions of multipolarity J'^ = 0^ — 5^. The in- 
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T[MeV] 

Fig. 7. Inclusive cross sections for the reaction e“)^®Co, averaged over the Fermi-Dirac dis¬ 

tribution, shown as a function of temperature in the interval T = 2 — 10 MeV. The RHB-I-RQRPA 
(DD-ME2) results are compared with those obtained using the shell model -|- RPA,I^ and the 
hybrid model of Ref. as well as those of RefS 


elusive cross sections and a full set of flux-averaged cross sections for the Z = 8 — 82 
pool of nuclei is available for the range of temperatures T = 0 — 10 MeV!^ As an il¬ 
lustrative example, the (r'e, e~) cross sections, averaged over the supernova neutrino 
flux ([3]) for T = 4 MeV, a = 0 are shown in Fig. [8] The cross sections are addition¬ 
ally marked for different groups of target nuclei: (i) stable nuclei, (ii) nuclei with 
proton excess N/Z < 1, and (iii) neutron-rich nuclei with with N/Z > 1.5. In the 
case of neutron-rich nuclei the cross sections are enhanced (blue squares) in com¬ 
parison to stable nuclei (up to an order of magnitude), because of the larger number 
of neutrons that can participate in charge-exchange neutrino-induced reactions. On 
the other side of the valley of stability, for neutron-deficient and proton-rich nu¬ 
clei, the {i'e,e~) cross sections are considerably reduced (crosses) due to the Pauli 
blocking of proton orbitals available for neutrino-induced transitions. 

The calculated neutrino-nucleus cross sections can be used for modeling the de¬ 
tector response to neutrinos involving different target nuclei, as well as in r-process 
calculations involving various astrophysical scenarios. It is, therefore, important to 
determine systematic uncertainties in the calculation of neutrino-nucleus cross sec¬ 
tions when various theoretical frameworks and effective interactions are employed. 
For the example of large-scale modeling of cross sections averaged over the super¬ 
nova neutrino flux for T=4, a = 0, the results of self-consistent calculations based 
on the RNEDF framework have been compared to those obtained using the RPA 
(WS-fLM)Pand ETFSI-fCQRPA.lSil Figure [9] displays the ratio of supernova neu- 
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Fig. 8. RHB+RQRPA (DD-ME2) inclusive neutrino-nucleus cross sections averaged over the 
Fermi-Dirac distribution for T = 4 MeV, a = 0, as a function of the mass number of target 
nuclei. The cross sections for particular groups of target nuclei are further marked with filled (red) 
circles for stable nuclei, crosses for nuclei with N/Z < 1, and filled (blue) squares for nuclei with 
N/Z > 1.5. 


trino flux-averaged cross sections calculated with the RHB-I-RQRPA (DD-ME2) 
and ETESI-|-CQRPA,l^as a function of the mass number of target nuclei. To iden¬ 
tify differences specific to particular groups of nuclei, the results are additionally 
labeled for stable nuclei, neutron-deficient nuclei with N/Z < 1, and neutron-rich 
nuclei with N/Z > 1.5. The RHB-I-RQRPA cross sections include higher multipoles 
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Fig. 9. The ratio of the RHB+RQRPA cross sections shown in Fig. |8] and the cross sections 
calculated using the ETFSI+CQRPApS framework. 


up to J = 5, whereas the ETFSI+CQRPA results contain only IAS and GT transi¬ 
tions. For most nuclei the RHB-f RQRPA approach based on the RNEDF framework 
predicts systematically larger cross sections, up to a factor of « 4 — 6. In the case of 
neutron rich-nuclei {N/Z > 1.5), for which the absolute values of the cross sections 
are significant, the differences between the two models behave in rather system¬ 
atic way and the ratio of the cross sections increases with an almost linear mass 
dependence. In the mass region A = 50 — 100, for the set of N/Z < 1 nuclei the 
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RHB+RQRPA cross sections are smaller than the corresponding ETFSI+CQRPA 
ones. The reason is to a large extent the anomalous enhancement of the cross 
sections for neutron-deficient nuclei predicted by the ETFSI-I-CQRPA model (cf. 
Ref. l30l for more details). 

7. Beta-decay half-lives 

An important application of elementary charge-exchange transition studies is the 
evaluation of /3-decay rates of neutron-rich nuclei. In the /3“ decay process 

z^N + e + P (4) 

a neutron is transformed into a proton, and this additional proton enhances the 
Coulomb repulsion pushing the isobaric analogue state above the Q/j energy win¬ 
dow. Fermi transitions, therefore, do not contribute to the decay of neutron-rich 
nuclei, and it is the low-energy part of the Gamow-Teller spectrum that determines 
the decay rat^^ 

rWo,i _ 

W\JW’^ - 1 (Wo.i - Wf F{Z, W)C{W)dW. (5) 

Here, W is the emitted electron energy in units of the electron mass, Wq.i is the 
maximum electron energy, F{Z,W) is the Fermi function that corrects for the 
Coulomb field of the nucleus. C{W) is the shape factor which, when only allowed 
transitions are taken into account, depends only on the Gamow-Teller transition 
strength. In neutron-rich nuclei, however, high-lying neutron states are partially oc¬ 
cupied enabling first-forbidden, parity-changing transitions, such as the spin-dipole, 
to provide a significant contribution to the total decay rate. In that case the shape 
factor becomes a complicated combination of transition matrix elements and terms 
that include the energy dependence. 

As the contribution of the phase space available to leptons depends on the 
energy of the transition, for the evaluation of decay rates it is important to achieve 
a good description of the low-energy part of the spectrum. Thus, the T = 0 pairing 
interaction Eq. o is an essential ingredient for any self-consistent mean-field model 
that describes the decay properties of neutron-rich nuclei.!^ In several previous 
studies it was found that a single value of the pairing strength cannot be used for 
all isotopic chains or in different mass regions.^^HIHllini jjowever, with the assumption 
that the strength of the T = 0 pairing depends on the difference between the number 
of neutrons and proton^^ 

Vq =Vl + 2 ^a+b(N-Z) ’ 

it is possible to obtain a very good agreement with data over the entire nuclear 
chart. The results presented in this section have been obtained using the values: 
Vl = 160.0 MeV, Vd = 15.0 MeV, a = 7.2 and b = —0.3, adjusted to obtain the 
best possible description of available data on /3-decay half-lifes 
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neutron number 

Fig. 10. /3-decay half-lives of neutron-rich isotopes of silver (top) and cadmium (bottom), obtained 
with two theoretical approaches and compared to experimental values. The FRDM-I-QRPA values 
are denoted by (red) squares, RHB-hRQRPA (D3C*) results by (blue) diamonds. For most nuclei 
the error bars on data points are too small to be discernible. 


In Fig. [TO] we display /3-decay half-lives calculated with the finite-range droplet 
model (FRDM) + QRPAP^(red squares) and the microscopic RHB-I-RQRPA model 
(effective interaction D3C*)P^ (blue diamonds), in comparison to available data.^ 
In the top panel results for the isotopes of silver {Z = 47) are shown for N = 67—90. 
A prominent feature is the relatively large discrepancy between model predictions 
for nuclei close to stability. The RHB-I-RQRPA (D3C*) results reproduce the em¬ 
pirical half-lives and follow closely the decreasing trend of the data, except in the 
case of ^^^Ag. The half-life of this nucleus is shorter by two orders of magnitude 
compared to neighboring isotopes, and this is very difficult to reproduce in model 
calculations. The FRDM reproduces the decay rate of ^^'^Ag but systematically un¬ 
derestimates the decay half-lives of longer-lived isotopes. At the very neutron-rich 
side both models predict results in good agreement with data. We notice, however, 
that the FRDM results display a sudden deviation from the general trend (between 
N = 86 and N = 90) by as much as an order of magnitude, while the D3C*-based 
RHB-fRQRPA predicts a very smooth mass dependence of the decay half-lives. 

In the case of the cadmium isotopic chain {Z = 48), the predictions of the two 
models are in better agreement and, generally, both models reproduce the data. 
Systematic deviations from the experimental values and RHB-I-RQRPA results are 
apparent in the FRDM predictions for the isotopes N > 76. One may notice that for 
both silver and cadmium isotopes the FRDM results display an odd-even staggering, 
not observed in the data. This staggering and the systematic deviation from the 
data may have important consequences in r-process simulations and the resulting 
elemental abundances. 
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In both isotopic chains the main contribution to the decay rate originates from 
the very strong Gamow-Teller transition 1 ^ 57/2 —>■ 7159/21 especially for nuclei closer 
to stability. With the addition of more neutrons, in particular above the N = 82 
closed neutron shell, additional configurations become available and first-forbidden 
transitions produce a significant contribution to the total decay rate. If Fig [Til we 



Fig. 11. Contribution (in percentage) of first-forbidden transitions to the total /3-decay rates of 
silver and cadmium isotopes. 


display the contribution of first-forbidden transitions to the total decay rate of silver 
and cadmium isotopes, from fV = 66 to = 92. Below the N = 82 closed neutron 
shell the contribution of first-forbidden transitions is rather small (up to 10%), 
and limited to the = 1“ vhxxji —>■ 7159/2 configuration. Above the closed shell, 
however, the neutron orbital / 7/2 and, to a lesser extent / 19 / 2 , become partially 
occupied and contribute to both the 0“ and the 1“ transitions. The occupation of 
opposite parity states in the next neutron shell leads to the rapid increase of the 
contributions shown in Fig. [TT] above the neutron shell closure. 

An important property of ,5-unstable nuclei is the number of neutrons emitted 
after the decay, that is, the /3-delayed neutron emission probability. The difference 
between the time scales of the /3-decay and the subsequent particle emission pro¬ 
cess justifies the assumption of their statistical independence. Thus, the /3-delayed 
neutron emission is considered as a multi-step process that proceeds starting with 
the /3-decay of the parent nucleus (A, Z), followed by either the 7 de-excitation of 
the daughter nucleus (A, Z 1), or by the emission of a neutron and the formation 
of the final nucleus {A—1,Z + 1). The probability of emission of a single or multiple 
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neutrons can be expressed 

^tot 

with the assumption that the daughter nucleus will emit as many neutrons as 
energetically allowed. Here, Xi is the partial decay rate, Xtot the total decay rate, 
and the sum takes into account all transitions with energies between x and (a; + 1) 
neutron separation energies. 




Fig. 12. Comparison of beta-delayed neutron emission probabilities calculated in the FRDM -|- 
QRPA (red squares) and RHB-hRQRPA (D3C*) (blue diamonds) models. The top panels display 
the one-neutron emission probabilities Pin for the silver (left) and cadmium (right) isotopic chains, 
while in the bottom panels the two-neutron emission probabilities P 2 n are shown. 


In Fig.[T2]we compare the one (Pin) and two (P 2 n) neutron emission probabil¬ 
ities for silver and cadmium isotopes. The odd-even staggering of the probabilities 
is a direct consequence of the staggering of the neutron separation energies in the 
daughter nuclei. A nucleus with an even number of neutrons decays into a nu¬ 
cleus with an odd number of neutrons, i.e. the one-neutron separation energy Sn 
is smaller in the daughter nucleus. In the case of a parent with an odd neutron 
number, the decay is into a daughter nucleus with an even number of neutrons and 
a correspondingly larger Sn- For a parent with an even number of neutrons the sum 
in Eq. Q takes into account a wider energy interval and this results in a higher 
neutron emission probability. 

For the lighter isotopes of both chains, the two models (FRDM -|- QRPA (red 
squares) and RHB-I-RQRPA (D3C*) (blue diamonds)) predict a similar dominance 
of the one-neutron emission up to V = 84. Significant differences are found above 
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N = 84 in silver isotopes. The FRDM predicts very small values of Pin, and 
comparatively much larger values of P 2 n for nuclei up to = 92, whereas the 
RHB+RQRPA (D3C*) yields a more gradual decrease of both Pin and P 2 n- Above 
N = 92 both models give small values of Pin and P 2 n, because for very neutron- 
rich nuclei the emission of 3 or more neutrons dominates. The differences between 
model predictions are much less pronounced in the case of Cd isotopes, especially 
for the two-neutron emission probability. 


8. Conclusions 

Weak-interaction processes provide valuable information on ground-state properties 
and the structure of excited states of a large number of nuclei that lie on both sides 
of the valley of ,5-stability. Calculations of stellar nucleosynthesis, nuclear aspects 
of supernova collapse and explosion, and neutrino-induced reactions, necessitate 
as input properties of thousands of nuclei far from stability, including characteris¬ 
tics related to weak-interaction processes. Many of these nuclei, especially on the 
neutron-rich side, are not accessible in experiments and their properties, therefore, 
must be determined using different theoretical approaches. One of these is the self- 
consistent mean-field method based on universal nuclear energy density functionals 
or global effective interactions, which can be applied in the description of arbitrar¬ 
ily heavy nuclei, exotic nuclei far from stability, and even nuclei at the nucleon 
drip-lines. 

We have reviewed some recent applications of a theoretical framework based 
on relativistic energy density functionals in modeling neutrino-nucleus reactions 
and calculation of /3-decay rates. Starting from a self-consistent mean-field solution 
for the nuclear ground state, a consistent proton-neutron relativistic quasiparti¬ 
cle random-phase approximation is used to calculate collective excited states in 
daughter nuclei and compute cross sections and transition rates. The model has 
very successfully been applied not only to the weak-interaction processes described 
in the present article, but also to muon capture and stellar electron capture.^^ 
The principal advantage of using this approach is the ability to perform system¬ 
atic calculations of weak-interaction rates in nuclei of arbitrary mass along the r- 
process path, based on microscopic self-consistent mean-field potentials rather than 
empirical potentials, and the inclusion of transitions of arbitrary multipolarity. Fu¬ 
ture developments will include the exploration of more complex configurations, e.g. 
particle-vibration coupling effects, a systematic treatment of finite-temperature ef¬ 
fects, and the extension of the RHB-I-RQRPA model to explicitly take into account 
quadrupole deformations. 
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